% devuelve un histograma de activacion promedio de 'tries' pruebas para 'M' inputs

function h = histograma(M,eta,tries)
	B = ones([M tries]);
	for i=[1:tries]
		%datos_entrenamiento = rand(M,1);
		datos_entrenamiento = normrnd(.5,.0625,M,1);
		datos_testeo = normrnd(.5,.0625,M,1);
		%datos_testeo = rand(M,1);
		
		red = som_lininit(datos_entrenamiento,'rect','sheet','msize',[M 1]);
		[red estado] = som_seqtrain(red,datos_entrenamiento,'trainlen',50,'alpha_ini',eta,'tracking',3);
		B(:,i) = som_bmus(red,datos_testeo);
	end
	H = hist(B,M);
	h = zeros(M,1);
	for i=[1:tries]
		h(:,1) = h(:,1) + H(:,i);
	end
	h = h/tries;
end
